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INTRODUCTION 

The  current  program  is  a  continuation  of  our  previous  studies  on  the  “computer  engineer¬ 
ing”  of  ionic  molecular  solids.  The  basic  aim  of  the  program  is  the  predictions  of  physical 
properties  by  First  Principles  calculations  entirely  free  from  disposable  parameters. 

Our  principal  theoretical  tool  is  the  Gordon-Kim  form  of  the  modified  electron  gas 
approximation,  supplemented  by  quantum  chemistry  for  the  treatment  of  covalently  bonded 
molecular  groups,  e.g.,  SOI”.  Most  recently  this  has  been  supplemented  by  the  self-consistent 
atomic  deformation  model  (SCAD)  and  full  augmented  plane  wave  band  structure  calcula¬ 
tions.  In  addition,  studies  have  been  made  on  the  phase  diagrams  of  binary  mixtures  of 
different  fluorides.  One  material  we  have  particularly  emphasized  is  sodium  calcium  triflu¬ 
oride.  While  this  has  only  been  synthesized  in  the  computer,  its  predicted  properties  are 
so  interesting  that  actual  physical  synthesis  could  produce  a  material  with  extremely  useful 
properties.  Specifically,  we  find  it  to  assume  the  lithium  niobate  structure  with  a  ferroelec¬ 
tric  polarization  at  ambient  temperature  comparable  to  barium  titanate.  This  suggests  that 
it  is  a  candidate  for  ferroelectric  applications  which  could  be  superior  to  conventional  oxide- 
based  systems.  In  particular,  it  is  transparent  well  into  the  ultra  violet  (UV)  and  could, 
for  example,  provide  a  non-linear  material  for  UV  optics,  as  distinct  from  oxides  which  are 
opaque  in  this  spectral  range. 

Since  the  abstracts  of  our  relevant  publications  contain  a  full  summary  of  each  inves¬ 
tigation,  we  shall  use  these  interspersed  with  appropriate  introductory  material. 

Studies  on  Mixed  Fluorides 

These  are  described  in  the  following  two  abstracts.  In  one  case,  KCaF3  is  mixed  with 
RbCaF3,  in  the  other  the  components  are  NaMgF3  and  KMgF3. 

As  can  be  seen  agreement  between  theory  (Gordon-Kim)  and  experiment  is  good,  in 
particular  an  observed  eutectic  point  in  the  KMgF3,  NaMgF3,  phase  diagram  is  reproduced 
by  theory. 

Structural  Phase  Transitions  in  Ionic  Molecular  Solids 

This  is  a  general  review  paper  which  discusses  at  length  the  basic  principles  which 
govern  the  origins  of  incommensurate  phases  in  these  systems.  It  brings  together  a  wide 
range  of  our  earlier  results  in  order  to  exemplify  these  principles. 

Structural  Phase  Transitions  in  Systems  with  Linear  Molecular  Ions 

These  papers  discuss  the  order-disorder  transitions  in  alkali  azides,  thiocyanates,  cyanides 
and  thallium  azide.  All  of  these  compounds  are  either  distorted  NaCl  or  CsCl  structures. 
The  transitions  involve  hindered  rotations  of  the  linear  ions  producing  a  variety  of  high  tem¬ 
perature  phases,  some  of  which  while  predicted  theoretically,  are  pre-empted  by  melting  in 
the  real  systems. 
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Studies  on  Mixed  Fluorides  Abstract 
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PHASE  TRANSITIONS  IN 

MIXED  ALKALI  CALCIUM  TRIFLUORIDE  SOLID  SOLUTIONS 
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(Refereed) 
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ABSTRACT 

Both  KCaFj  and  RbCaFj  are  known  to  undergo  ferroelastic  phase  transitions, 
the  former  at  temperatures  near  550  K  and  the  latter  at  195  and  40  K.  The  two 
compounds  comprise  ends  of  the  completely  miscible  solid  solution 
Ki-sRbxCaFj.  Simulations  of  molecular  dynamics  have  been  performed  for 
this  solid  solution,  where  Ro*  ions  were  substituted  for  K*  ions  in  KCaFj  in 
varying  concentrations  of  randomly  distributed  defects  representing 
constitutional  crossover  to  RbCaFj.  First-principles  calculations  were  used  to 
obtain  the  potentials  between  ion  pairs  in  these  samples;  the  pair  potentials 
were  used  to  calculate  the  temperature  dependence  of  the  interaxis  angles 
which  were  used  as  markers  for  the  polymorphic  phase  transition.  The 
calculated  transition  temperature  thus  obtained  decreases  in  a  nearly  linear 
fashion  as  a  function  of  Rb+  ion  concentration.  Experimental  transition 
temperatures  obtained  from  differential  scanning  calorimetry  of  actual 
samples  agree  well  with  the  calculated  temperatures  for  those  compositions 
with  transitions  above  room  temperature. 


KEYWORDS:  A.  fluorides,  C.  differential  scanning  calorimetry’  (DSC),  D. 
phase  transitions 
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Abstract 

Phase  transitions  in  the  miscible  solid  solution  Nal-xKxMgF5  were  examined  over  a  wide  range 
of  compositions  by  computer  molecular  dynamics.  X-ray  diffraction,  and  thermal  analysis  in  order 
to  characterize  the  polymorphic  phase  transitions  as  a  function  of  alkali-metal  content.  The  pure 
NaMgF5  composition  has  a  single  orthorhombic-to-cubic  transition  at  1038  K.  but  computer 
modeled  compositions  with  K~  partially  substituted  for  Na*  ions  have  at  least  two  polymorphic 
transitions.  The  models  also  indicate  that  the  transition  temperatures  decrease  with  increasing 
potassium  content.  Results  from  thermal  analyses  and  from  literature  give  similar  results. 
Computer  simulations,  experimental  data,  and  literature  values  all  show  a  room-temperature 
transition  for  the  composition  around  NaQ  65Kq  35MgF5.  ©  2000  Elsevier  Science  Ltd.  All  rights 
reserved. 

Keywords:  A.  Fluorides:  C.  Differential  scanning  calorimetry  (DSC):  C.  X-ray  diffraction:  D.  Lattice  dynamics: 
D.  Phase  transitions 


*  Corresponding  author.  Tel.:  -r 402-554-3592:  fax:  -M02-554-388S. 
E-mail  address:  robert.smith@unomaha.edu  (R.W.  Smith). 

O025-54f)3/{X)/S  -  see  front  matter  €  2000  Elsevier  Science  Ltd.  All  rights  reserved. 
Pli:  S002 5 -5408(00)0023 2-4 
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STRUCTURAL  PHASE  TRANSITIONS  IN 
IONIC  MOLECULAR  SOLIDS* 

J.R.  HARDY 

Department  of  Physics  and  Center  for  Electro-Optics.  University 
of  Nebraska-Lincoln .  Lincoln .  NE  68588-0111,  USA 


An  overview  is  presented  of  our  studies  on  the  nature  of  structural  instabilities  in  rela¬ 
tively  complex  ionic  solids.  These  are  based  on  parameter-free  interionic  potentials  based 
on  the  Gordon- Kim  modified  electron  gas  formalism  extended  to  molecular  ions. 

We  describe  the  manner  in  which  there  emerge  from  these  studies  quite  general  con¬ 
cepts  of  “size”  and  “shape"  as  structural  determinants.  In  particular,  we  discuss  how 
these,  and  the  approximate  symmetries  that  they  can  produce,  can  provide  a  relatively 
simple  structure- based  explanation  of  the  origins  of  incommensurate  phases  in  these 
systems.  However,  we  also  emphasize  that  the  existence  of  such  symmetries  does  not 
guarantee  an  incommensurate  phase.  This  can  only  be  realized  if  long-range  correlations 
are  sufficiently  strong  to  overcome  random  local  disordering.  Thus,  either" the  molecular 
units  are  partially  linked  and/or  there  exist  long-range  Coulomb  interactions  between 
individual  units. 

Keywords:  Lattice  dynamics:  Ionic  molecular  solids;  Incommensurate  structures 
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Perhaps  the  most  interesting  system  is  KSCN.  Experimentally,  this  appears  to  assume 
a  tetragonal  structure  just  prior  to  melting.  However,  mircoscopically,  it  is  found  to  consist 
of  dynamically  moving  domains  of  the  low-temperature  orthorhombic  phase.  Theoretically, 
our  molecular  dynamics  supercell  is  far  too  small  to  reveal  this  behavior.  However,  what  we 
do  observe  is  the  onset  of  CNS-  flipping  out  of  the  basal  plane.  We  argue  that  this  is  the 
origin  of  the  moving  domain  walls  producing  the  dynamically  twinned  structure  observed 
experimentally. 

One  of  other  significant  finding  is  that  our  combined  quantum  chemistry  Gordon-Kim 
approach  produces  activation  barriers  for  CN“  rotations  that  are  far  too  high,  indicating 
that  a  more  sophisticated  approach  (e.g.,  SCAD)  is  necessary  for  these  strongly  polar  ions. 

More  Sophisticated  First  Principles  Studies 

The  first  two  investigations  under  this  heading  are  applications  of  the  self-consistent 
atomic  deformation  model  (SCAD).  In  this  model,  localized  atomic  states  centered  on  the 
respective  ions  are  used.  Each  is  then  calculated  self-consistently  in  the  crystal  field  of  the 
other  ions.  Thus,  as  the  lattice  is  deformed,  the  ions  deform  in  a  self-consistent  manner; 
the  deformation  at  each  site  being  consistent  with  the  field  at  that  site  due  to  the  other 
(deformed)  ions.  Hence,  the  overall  charge  distribution  and  total  energy  change  are  self- 
consistent. 

In  the  first  study,  the  method  is  applied  to  the  alkali  halide  sequence  with  very  rea- 

\ 

sonable  results:  in  particular,  polar  motions  produce  Bom  effective  charges  in  good  accord 
with  experiment  as  distinct  from  the  simple  Gordon-Kim  result  of  unity  in  all  cases. 

In  the  second  study,  the  method  is  applied  to  NaCaF3  to  examine  the  effect  of  this 
modification  on  the  ferroelectric  behavior  of  this  system.  Once  again,  the  system  is  found  to 
be  ferroelectric  with  a  well  depth  comparable  to  the  Gordon-Kim  value.  Also,  the  compound 
remains  metastable  with  respect  to  the  constituents  NaF  and  CaF2. 

Finally,  we  present  a  first  principles  study  of  the  optical  (linear  and  non-linear)  prop¬ 
erties  of  KNbOs. 

Further  Studies  on  Ionic  Molecular  Solids  within  the  Gordon-Kim  Lu-Hardy 
Approach 

This  approach  is  that  employed  in  our  previous  work  on  ionic  molecular  systems.  The 
molecular  ion  charge  density  is  calculated  using  the  GAUSSIAN  Quantum  Chemistry  package 
and  it  is  then  decomposed  into  individual  atomic  components  which  are  employed  in  a 
Gordon-Kim  derivation  of  all  potentials  external  to  the  group.  The  internal  motion  is  treated 
within  the  harmonic  approximation  using  force  constants  generated  by  GAUSSIAN.  For 
the  monatomic  cations,  charge  densities  calculated  from  the  tabulated  wave  functions  are 
employed  in  the  respective  Gordon-Kim  calculations. 
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PHYSICAL  REVIEW  B  VOLUME  60,  NUMBER  22  1  DECEMBER  1999-11 

Molecular-dynamics  study  of  phase  transitions  in  alkali  azides 

M.  M.  Ossowski  and  J.  R.  Hardy 

Department  of  Physics  and  Center  for  Electro-Optics .  University  of  Nebraska.  Lincoln.  Nebraska  68588-02 JJ 

R.  W.  Smith 

Department  of  Chemistry.  University  of  Nebraska.  Omaha.  Nebraska  68182-0109 
(Received  18  June  1999) 

An  account  is  presented  of  our  studies  of  the  order-disorder  phase  transitions  in  KN3,  RbX>  and  CsN3. 

These  are  based  on  parameter-free  interionic  potentials  based  on  the  Gordon-Kim  modified  electron-gas  for¬ 
malism  extended  to  molecular  ions.  We  performed  static  structural  relaxations  and  supercell  molecular  dynam¬ 
ics  and  predicted  with  reasonable  accuracy  the  temperatures  for  the  onset  of  the  transitions.  In  particular,  we 
address  the  question  of  how  the  N3~  ions  reorient  to  yield  the  transitions.  We  found  the  existence  of  NaCl-type 
high-temperature  phases  in  disordered  KN3  and  RbN3  and  argue  that  this  restructuring  is  preempted  by  melting 
in  these  two  systems.  [SO  163 -1829(99)  1054 1-1] 


PHYSICAL  REVIEW  B  VOLUME  62,  NUMBER  5  1  AUGUST  2000-1 

Molecular-dynamics  study  of  phase  transitions  in  alkali  thiocyanates 

M.  M.  Ossowski  and  J.  R.  Hardy 

Department  of  Physics  and  Center  for  Electro-Optics .  University  of  Nebraska,  Lincoln .  Nebraska  68588-0111 

R.  W.  Smith 

Department  of  Chemistry.  University  of  Nebraska.  Omaha .  Nebraska  68182-0109 
(Received  22  November  1999) 

An  account  is  presented  of  our  studies  of  the  .order-disorder  phase  transitions  in  KSCN.  RbSCN.  and 
CsSCN.  These  are  based  on  parameter-free  interionic  potentials  based  on  the  Gordon-Kim  modified  electron 
gas  formalism  extended  to  molecular  ions.  We  performed  static  structural  relocations  and  supercell  molecular 
dynamics  and  predicted  with  reasonable  accuracy  the  temperatures  for  the  onset  of  the  transitions.  In  particular, 
w'e  address  the  question  of  how  the  SCN~  ions  disorder  to  yield  subsequent  structural  transformations.  We 
found  high-temperature  phases  of  average  Fm?m  symmetry  for  both  KSCN  and  RbSCN.  We  argue  that  in 
reality  the  full  appearance  of  these  phases  is  preempted  by  melting.  However,  they  are  candidates  for  the  twin 
boundaries  which  are  observed  in  the  high-temperature  “average"  tetragonal  phases.  The  high-temperature 
phase  of  CsSCN  was  found  to  be  of  average  Fmim  symmetry'. 
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Computer  Molecular  Dynamics  and  Phase  Transitions  in  Alkali  Azides  and 

Thiocyanates 

M.  M.  Ossowski  and  J.  R.  Hardy 

Department  of  Physics  and  Center  for  Electro-Optics,  University  of  Nebraska,  Lincoln,  Nebraska  68588- 

0111,  USA,  martin@solids.unl.edu 


ABSTRACT 

An  account  is  presented  of  our  studies  of  the  order- 
disorder  phase  transitions  in  KNj,  RbNj,  CsNj,  KSCN, 
RbSCN  and  CsSCN.  These  are  based  on  parameter-free 
inter-ionic  potentials  based  on  the  Gordon-Kim  modified 
electron  gas  formalism  extended  to  molecular  ions. 

With  these  potentials  \ve  performed  static  structural  re¬ 
laxations  and  supercell  molecular  dynamics  and  predict 
with  reasonable  accuracy  the  temperatures  for  the  onset  of 
the  transitions.  In  particular  we  address  the  question  of 
how  the  Nf  and  SChT  ions  reorient  to  yield  the  transitions. 

It  was  found  that  in  case  of  azides  we  observed  hin¬ 
dered  rotations  of  the  anions  about  all  three  crystallographic 
axes,  in  all  three  systems.  However,  in  the  thiocyanates 
only  CsSCN  yielded  hindered  rotations  about  all  three  axes. 
In  KSCN  and  RbSCN  the  order-disorder  transition  appears 
rather  to  involve  large  amplitude  librations  of  the  SCN" 
ions,  primarily  about  the  c  axis. 

Keywords:  Azides;  Thiocyanates;  Ionic  molecular  solids; 
Phase  transitions 
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Abstract 

The  phase  transition  in  T1N3  is  simulated  based  on  the  potentials  calculated  from  the  Gordon-Kim 
modified  electron  gas  model  extended  to  ion  molecular  crystals.  It  is  found  that  TIN3  transforms  into 
a  cubic  CsCl  structure  at  high  temperature  due  to  the  rotations  of  the  N3  ions.  Above  the  phase 
transition  the  orientations  of  the  N3“  ions  are  random  with  four  preferred  orientations  with  respect  to 
the  cubic  axes.  ©  2001  Elsevier  Science  Ltd.  All  rights  reserved. 

Keywords:  D.  phase  transitions 


Structural  Transitions  in  NaCN  and 

KCN 

M.  M.  Ossowski  ancl  J.  R.  Hardy*  and  R.  W.  Smith1 
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University  of  Nebraska.  Lincoln.  Nebraska  dS5SS-0Ul 
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Abstract.  An  account  is  presented  of  our  studies  of  the  order-disorder  pluisc  tran¬ 
sitions  in  XaCX  and  KGX.  Those  are  based  on  parameter-free  inter-ionic  potentials 
based  on  the  Gordon-Kiui  modified  electron  gas  formalism  extended  to  molecular  ions. 

We  performed  static  structural  relaxations  and  suporceli  molecular  dynamic's  juhI 
repnxluced  two  transitions  known  in  each  of  those  systems. 

We  also  calculated  upper  bouucls  to  the  barrier  of  rotation  of  a  cyanide  ion  in  a 
ground  state  of  XaCX  and  KCX  ami  discuss  possible  short  comings  of  our  model. 
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More  Sophisticated  First  Principles  Studies  Abstract 


PHYSICAL  REVIEW  B  VOLUME  6 ! .  NUMBER  1 7 


1  MAY  2000-1 


Calculation  of  electronic,  structural,  and  vibrational  properties  in  alkali  halides 
using  a  density-functional  method  with  localized  densities 

W.N.  Mei 

Department  of  Physics.  University  of  Nebraska  at  Omaha.  Omaha.  Nebraska  68182-0266 


L.  L.  Boyer  and  M.  J.  Mehl 

Center  for  Computational  Materials  Science.  Naval  Research  Laboratory.  Washington.  D.C.  20375-5245 

M.  M.  Ossowski 

Department  of  Physics.  University  of  Nebraska  at  Lincoln.  Lincoln.  Nebraska  685SS 


H.T.  Stokes 

Department  of  Physics  and  Astronomy.  Brigham  Young  University.  Provo.  Utah  84602 
(Received  7  December  1999) 

A  recently  developed  density-functional  method  based  on  localized  densities  is  applied  to  calculate  elec¬ 
tronic.  structural,  and  vibrational  properties  of  20  alkali  halides  with  elements  lithium  through  cesium  and. 
fluorine  through  iodine.  Properties  calculated  include  dissociation  energy,  lattice  parameter,  dielectric  constant, 
elastic  moduli,  and  phonon  frequencies  for  the  high-symmetry  points  of  the  Brillouin  zone.  Results  are  dis¬ 
cussed  and  compared  with  experiment  and  other  calculations. 
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More  Sophosticated  First  Principles  Studies  Abstract 


Predicted  Properties  of  NaCaFs 

L.  L.  Boyer  and  M.  J.  Mehl 

Center  for  Computational  Materials  Scicva:,  Kara!  Research  Laboratory.  Washington,  D.C. 

20375-5X45 

W.  N.  Mei.  Chun-gang  Duan.  J.  W.  Flocken  and  R.  A.  Guenther 

DciHirtmr.nl  of  Physics.  University  of  Nebraska  at  Omaha.  Omaha.  Nebraska  GSIS2-02GG 

J.  R.  Hardy 

Department  of  Physics  and  Center  for  Electro-  Optics.  University  of  Nebraska  at  Lincoln. 

Lincoln.  Nebraska  SS5SS 

H.  T.  Stokes 

Department  of  Physics  and  Astronomy.  Brigham  Young  University.  Provo.  Utah  #4602 

P.  J.  Edwardson 

Mission  Research  Corporalwiu  Nashua.  NH  03062 

Abstract.  More  than  a  <I<vade  ago  computer  simulations  based  on  Gordoti-Kiut  po¬ 
tentials  suggested  that  NaCaF*  would  be  a  ferroelectric  compound  if  it  could  be  formed 
experimentally.  Recent  attempts  to  form  thin  films  of  NaCaF*  using  pukod-laser  de¬ 
position  haw  prompted  us  to  carry  out  further  theoretical  studies  of  this  material. 
Here  \vc  apply  the  self-consistent  atomic  deformation  method  to  calculate  electronic, 
structural  and  vibrational  properties  of  NaCaFs  and  the  constituent  compounds  NaF 
and  CaFn. 
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More  Sophosticated  First  Principles  Studies  Abstract 
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First-principles  study  on  the  optical  properties  of 
KNb03 


Chun-gang  Duan\  W  N  Mci1,  Jianjun  Liu2  and  J  R  Hardy2 

1  Department  of  Physics.  University  of  Nebraska  at  Omaha.  Omaha.  NE  6SI  $2-0266.  USA 

2  Department  of  Physics  and  Center  for  Electro-Optics.  Universitv  of  Nebraska  at  Lincoln. 
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Online  at  stacks.iop.org/JPhysCM/13/S189 

Abstract 

We  report  our  studies  on  the  electronic  structure  and  linear  and  nonlinear  optical 
(NLO)  properties  of  KNbOy  using  a  first-principles  method  in  the  local  density- 
approximation  (LDA).  The  calculated  results  for  the  refractive  indices  and 
second-harmonic-generation  (SHG)  coefficients  agree  well  with  experimental 
results.  From  decomposing  the  nonlinear  susceptibility,  we  find  that  the  primary- 
contribution  to  the  NLO  behaviour  comes  from  the  hybridization  of  the  O  2p 
and  Nb  4d  electron  states.  In  addition,  there  are  two  different  roles  played  by 
the  0  atoms  because  of  their  different  distances  from  the  Nb  atom:  thus  we 
propose  a  possible  way  to  enhance  the  SHG  coefficients. 
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The  following  compounds  have  been  studied:  alkali  nitrites;  isomorphous  NaNC>3  and 
CaCC>3;  RbN03  and  CSNO3;  AgN03;  alkali  perdiolorates;  and  alkaline  earth  silicates. 

As  before,  this  work  is  detailed  in  the  following  abstracts  of  the  respective  studies. 
The  nitrites  are  of  particular  interest  insofar  as  the  NOJ  ion  has  a  permanent  dipole 
and  thus,  in  the  high  temperature  phases  where  these  ions  axe  in  quasi-free  rotation,  they 
can  couple  to  an  applied  microwave  field  and  provide  loss  and  optical  dispersion.  These 
compounds  thus  have  potential  applications  for  microwave  optics.  Of  the  other  compounds, 
the  silicates,  in  particular,  are  of  geophysical  interest  -  specifically  their  phase  diagrams  are 
the  subject  of  controversy  and  our  results  may  serve  to  settle  this  dispute. 

PERSONNEL 

Principal  Investigator:  Dr.  John  R.  Hardy 

Graduate  Student /Postdoctoral  Research  Associate:  Martin  Ossowski 
Visiting  Scholar:  Jianjun  Liu 


13 


Further  Studies  on  Ionic  Molecular  Solids  within  the  Gordon-Kim 
Lu-Hardy  Approach  Abstract 

PHYSICAL  REVIEW  B,  VOLUME  63,  144105 

Order-disorder  phase  transitions  in  KN02,  CsN02,  and  T1iNt02  crystals: 

A  molecular  dynamics  study 

Chun-gang  Duan*  and  \V.  N.  Mei 

Department  of  Physics,  University  of  Nebraska -Omaha,  Nebraska  6S  IS 2-0266 

R.  W.  Smith 

Department  of  Chemistry,  University  of  Nebraska  -  Omali  a,  Nebraska  6S1S2-0109 

Jianjun  Liu,  M.  M.  Ossowski,  and  J.  R.  Hardy 

Department  of  Physics  and  Center  for  Electro-optics,  University  of  Nebraska-Lincoln.  Nebraska  6S5SS-0IJ  J 
{Received  5  January  2000;  revised  manuscript  received  14  November  2000;  published  16  March  2001) 

The  order-disorder  phase  transitions  of  KN02,  CsN02,  and  T!NO:  have  been  studied  using  parameter-free 
molecular  dynamics  simulation.  It  is  found  that  the  phase  transitions  in  nitrites  investigated  are  driven  by  the 
rotations  of  the  N02~  ions  about  different  axes  together  with  displacements  of  cations  and  anions.  We 
successfully  reproduce  the  high-temperature  phases  of  these  nitrites,  i.e.,  the  NaCl-like  structure  for  KN02  and 
CsCl-Iike  structure  for  Cs(Tl)N02.  Based  on  the  investigation  of  the  radial  distribution  function  of  the  cations 
and  anions,  we  explain  why  KN02  and  Cs(Tl)N02  form  quite  different  low-temperature  phases. 

DOL:  10. 1  l03/PhysRevB.63. 144 105  PACS  numbers):  64.60.Cn,  6l.43.Bn,  64.70.Pf 


Phvs  Chem  Minerals  (2001)  23:  586-590 
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Jianjun  Liu  •  C.-G.  Duan  •  M.  M.  Ossowski 
W.  N.  Mei  *  R.  \V.  Smith  •  J.  R.  Hardy 

Simulation  of  structural  phase  transition  in  NaN03  and  CaC03 

Ricsivid:  30  January  2001  /  Accepted:  II  May  2001 


Abstract  The  order-disorder  phase  transitions  in  NaNO; 
and  CaCOj  are  simulated  by  molecular  dynamics.  The 
simulations  are  based  on  the  potentials  calculated 
from  the  Gordon-Kim  modified  electron  gas  formalism 
extended  to  molecular  ions.  We  successfully  reproduced 
the  transition  temperature  Tf  and  the  abnormally  large 
c  axis  thermal  expansion  observed  in  experiment.  The 
phase  transitions  in  NaNO;  and  CaCO;  were  found  to 
be  initiated  by  =  60  and  =  1 80s  reorientation  of  the 
NO;  and  CO;  ions  about  the  c  axis.  The  orientations  of 
NO;  and  CO;  ions  are  continuous  with  six  preferred 
calcite-type  orientations  above  the  phase-transition 
temperature. 


Key  words  Molecular  dynamics  simulation  • 

Gordon-Kim  potentials  •  Phase  transition 
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Further  Studies  on  Ionic  Molecular  Solids  within  the  Gordon-Kim 
Lu-Hardy  Approach  Abstract 

Simulation  of  Structural  Transformation 
in  Aragonite  CaCC>3 

Jianjun  Liu1,  M.  M.  Ossowski,  and  J.  R.  Hardy 


Dc/Htrtoteiu  of  Physic*  tout  Center  far  Electro-Optic*. 
IhthcruiyvfKehrtfiku,  Lincoln.  Nel*ru\lu  6SSSS-UI  /  /.  USA. 


Chun-gung  Duan  and  W.N.  Mei 


Department  of PUyttcs.  University' of  Nebraska.  Omaha.  A  cbnrda  6S 132-0266.  USA. 


Abstract.  The  structural  transformation  in  aragonite  CaCOj  is  simulated  hy  molecular 
dynamics.  llic  simulations  arc  Iwscd  on  the  potentials  calculated  from  the  Ciurdon-Kim  modified 
electron  t*:is  formalism.  We  found  two  pluisc  transitions  in  aragonite  at  high  temperature. 


lTM5.  fyn-loircul  F*iv«r$  cf  Ferret  leaner  20tX>  Xtprt  Center  for  Fhxvcs  Hmfrr  IVUfcp-.  edited  by  R.  E.  Cohen 
C  'i\0  \nrVK  !r.u;:*.rte  »»t  (0 

;>.>s 


Journal  of  Solid  State  Chemistry  160,  222-229  (2001)  ^ 

doi:I0.1006/jssc.200l.9226.  available  online  at  http://www.idealibrary.com  on  IDE 


Molecular  Dynamics  Simulation  of  Structural  Phase  Transitions 

in  RbNOs  and  CsN03 

Jianjun  Liu,*-1  Chun-gang  Duanj  M.  M.  Ossowski  *  W.  N.  Meij 
R.  W.  Smith,!  and  J.  R.  Hardy* 

"  Department  of  Physics  end  Center  for  Electro-Optics.  University  of  Nebraska.  Lincoln.  Nebraska  63533-01 1 1:  f  Department  of Physics.  University  of Ne brus< 
Omaha ,  Nebraska  6S 132-0266:  and  $  Department  of  Chemistry.  University  of  Nebraska.  Omaha.  Nebraska  63132-0109 

Received  January  12.  2001;  in  revised  form  April  IS.  2001;  accepted  April  30.  2001 


Structural  phase  transitions  in  RbNOj  and  CsNOy  are  studied 
by  molecular  dynamics.  The  simulations  are  based  on  the  para- 
meter-free  potentials  calculated  from  the  Gordon-Kim  modified 
electron  gas  formalism,  extended  to  ionic  molecular  crystals. 
The  microscopic  mechanism  of  the  structural  phase  transitions  in 
RbNCV,  and  CsNOj  is  revealed.  It  is  found  that  the  phase  IV— III 
transition  in  RbNOy  and  the  phase  II— I  transition  in  CsNOj  are 
initiated  by  the  in-plane  and  out-of-plane  rotations  of  the  NOj 
ions,  and  the  phase  III— II— I  transitions  in  RbNOj  are  due  to 
dilation  along  a  trigonal  axis  of*  phase  III,  giving  phase  II 
a  rhombohedral  structure,  z  2001  A«d«mc  pr«s 
Key  Words :  molecular  dynamics;  Gordon-Kim  potential; 
nitrate;  phase  transition;  order-disorder. 
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Further  Studies  on  Ionic  Molecular  Solids  within  the  Gordon-Kim 
Lu-Hardy  Approach  Abstract 
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Molecular  dynamics  simulation  of  phase  transition  in  AgNOs 

Jianjun  Liu®  *,  Churi-gang  Duanb,  M.M.  Ossowski®,  W.N.  Meib,  R.W.  Smith®, 

J.R.  Hardy® 
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Abstract 

Structural  phase  transition  in  AgNOj  at  high  temperature  is  simulated  by  molecular  dynamics.* The  simulations  are  based  on 
the  potentials  calculated  from  the  Gordon-Kim  modified  electron-gas  formalism  extended  to  molecular  ionic  crystals.  AgNO; 
transforms  into  rhombohedral  structure  at  high  temperature  and  the  phase  transition  is  associated  with  the  rotations  of  the  NO.s 
ions  and  displacements  of  the  N03  and  Ag  ions.  ©  2002  Elsevier  Science  Ltd.  All  rights  reserved. 

Keywords:  A.  Inorganic  compounds:  C.  Ab  initio  calculations:  D.  Phase  transitions 
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Order-Disorder  Structural  Phase  Transitions  in  Alkali  Perchlorates 

Jianjun  Liu.*T-1  Chun-gang  Duan,t  W.  N.  Mei.i  R.  W.  Smith.;  and  J.  R.  Hardy* 
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Received  June  15.  2001:  in  revised  form  September  20.  2001:  accepted  September  27.  2001 


Order-disorder  structural  phase  transitions  in  alkali  perchlor¬ 
ates  MC\Oa  (M=  Na,  K,  Rb,  Cs)  are  investigated  using  molecu¬ 
lar  dynamics  simulation.  The  potentials  in  the  simulations  are 
based  on  the  Gordon-Kim  modified  electron  gas  formalism 
extended  to  molecular  ions.  The  simulations  yield  first-order 
phase  transitions  in  perchlorates  from  low  temperature  ortho¬ 
rhombic  structures  to  high  temperature  cubic  NaCl  structures. 
The  perchlorate  ions  are  found  to  be  orientational  disordered  in 
the  high  temperature  phases.  «  aw:  Science 

Key  Words:  alkali  perchlorate;  molecular-dynamics  simula¬ 
tion;  phase  transition;  Gordon-Kim  potential;  order-disorder. 
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Structural  phase  transitions  in  Ca2Si04  and  Sr2Si04  are  investigated  by  molecular  dynamics 
simulations.  The  simulations  are  based  on  the  potentials  calculated  from  the  Gordon-Kim  modified 
electron  gas  formalism  extended  to  molecular  ions.  We  successfully  reproduced  the  transition  y 
-  aH~  a  and  fi - a'L  - a  in  Ca2Si04,  and  the  transition  f-a'  in  Sr2Si04.  We  find  that  the  a'L 
phase  of  Ca2Si04  is  an  aX3bXc  superstructure  of  the  a'H  phase,  while  the  ct'H  phase  has  a 
£-K,S04  structure,  and  the  a  phase  of  Ca2Si04  has  a  disordered  structure  with  space  group 
P6-Jmmc.  ©  2002  American  Institute  of  Physics.  [DOI:  10.1063/1.1446043] 
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